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We report on the serendipitous discovery of two additional 
6.7 GHz methanol masers in the same region, not associated 
with IRS 1 . Interferometry maser positions are compared with 
recent single-dish and interferometry continuum observations. 
The positions of the masers agree to high accuracy with the 
1.2 mm continuum peak emission in NGC7538 IRS 9 and 
NGC7538 S. This clear association is also confirmed by the po- 
sitional agreement of the masers with existing high resolution 
continuum observations at cm and/or mm wavelengths. Making 
use of the established strong relation between methanol masers 
and high-mas star formation, we claim that we have accurately 
positioned the high-mass protostars within the regions where 
they are detected. The variety of objects hosting a 6.7 GHz 
methanol maser in NGC 7538 shows that this emission prob- 

' ably traces different evolutionary stages within the protostellar 

■ phase. 
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Abstract. NGC7538 is known to host a 6.7 and 12.2 GHz methanol maser cospatial with a Ultra Compact (UC)Hn region, 
IRS 1. 



1. Introduction 

The NGC 7538 nebula is a furnace of intense massive star for- 
mation located at some 2.7 kpc from the Sun. At least 1 1 in- 
frare d sources ha ve been identified embedded within it (Fig.^ 
lOiha et aTll2004l and references therein). These IR sources lie 
in four regions hosting a large range of different types of young 
stellar objects (YSOs) with signs of decreasing age moving 
from the North- West to the South-East. In the far North- West 
the first star formation region contains the IR sources IRS 5, 6, 
7 which are associated with a developed H n region o f ~ 3 pc in 
size. IRS 6 has been proposed as its exciting source (Oiha et all 
120041) . A further star formation region can be identified with 
the IR sources IRS 1-3, where infrared emission at K-b and 
and l ong-ward is dominated by IRS 1 foe Buizer & MinieJ 
120051) . IRS 1 is known to power a Ultra Compact (UC) Hn re- 
gion, an early stage of massive star formation. Masers of dif- 
ferent species have also been detected toward this source (OH, 
H2O, CH3OH , see iMinier et alJll998l and references therein). 
A third star-forming region is centered on NGC 7538 S and 
also contains IRS 11, 10" to its northwest (Fig.Q. Coincident 
with NGC 7538 S O H and H 2 masers have been detected 
jKamevaet alj fl990). Finally, furthest to the South-East, the 
fourth known region of star formation contains the deeply em- 
bedded IR source IRS 9, believed to be in a stage prior to the 
formation of a UC H 11 region. H?Q masers were also dete cted 
toward this source foavis et alJl998t lKameva et al. 1990). 

Methanol masers are indicative of the earliest stages of 
massive star formation. They arise from cold, embedded 
molecular clumps and in some cases from dark infrared clouds 
Jffill et alJl2005l IPurcell et al]l2005l) . These clumps are often 
interpreted to be protoclusters of high mass YSOs, i.e. pre- 
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curso rs of OB associations (e.g. lMotte et alJ2003UMinier et alJ 
l2005h . High angular resolution observations of methanol maser 
environments reveal that these masers are associated with hot 
cores, UC H 11 regions and very often do not coincide with ra- 
dio continuu m emission fr om bright Ult ra Co mpact (UC) H n 
regions (Wals h et aDll998l IMinier et alJl200ll) . The detection 
of a methanol maser in such regions is a clear indication of the 
presence of a deeply embedded hot core that is characterised by 
large methanol abundance and 1 00-200 K gas and dust temper- 
ature. These physical conditions are in agreement with t hose 
derived from maser modelling, as e.g. Sob olevet all 19971 The 
high astrometric accuracy of VLBI observations of methanol 
masers gives thus the best estimate of the position of the high- 
mass protostellar core at a scale of ~ 10 AU. One clear exam- 
ple of this fact is given by the main component of the CH3OH 
maser in NGC 7538 (component A in Fig. 0: the maser pin- 
points the massive protostellar core within the IRS 1 submm 
clump to mill iarcsecond accuracy (1 mas at ~3kpc is ~3 AU, 
see lPestalozzi et all2004l) . 

In this letter we report on the serendipitous detection of 
6.7 GHz methanol masers in the third and fourth star formation 
regions in NGC 7538, cospatial with NGC 7538 S and IRS 9. 
Supported by what exposed above, we argue that with this dis- 
covery we have accurately positioned the two massive proto- 
stars in those regions. This is then confirmed by the coincidence 
of the new masers with the peak of the 1.2 mm dust emission 
coming from the thick cocoon surrounding them. 

2. Observations and data reduction 

The data used for this paper was taken in two runs, with the 
principal aim of studying the known 6.7 GHz methanol maser 
emission toward NGC 7538 (i.e. spectral components A-E, see 
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Fig. EJ- As explained below, the results presented here refer to 
the serendipitous discovery of two new maser features, unseen 
until now in any interferometry experiment. 

The first data set (runl, November 4"' -5"' 2004) was ob- 
tained during a joint experiment including the European VLBI 
Network 1 and the three antennas of the MERLIN 2 array which 
were then operational at 5cm wavelength, Cambridge (CM), 
Knockin (KN) and Darnhall (DA). The spectral resolution was 
1.9 and 3.9 kHz/chan (or 0.087 and 0.175 kms 'per chan- 
nel, respectively) over 512 channels for the joint EVN and the 
MERLIN only data, respectively. Only the MERLIN portion of 
this data set has so far been reduced and is presented here. The 
second data set (run2, December 11'*- 12* 2004) was a 13 hour 
MERLIN only experiment, including 5 antennas: CM, KN, 
DA, the 25 -m Jodrell Bank Mk2 antenna (JB) and Tabley (TA). 
The spectral resolution was 0.976 kHz/chan (0.044 kms~'per 
channel) over 512 channels. In both runs the telescopes were 
pointed, and the data correlated, at the position of the brightest 
methanol maser feature in IRS 1 , feature A in Figure[2](see also 
Table!]}. 

First order calibration and translation into fits format was 
performed with tdproc at Jodrell Bank. Corrections to the pri- 
mary amplitude calibration, phase calibration (self calibration 
on the brightest emission channel) and image production were 
performed in AIPS. To locate all the maser emission, large 
10x10 arcsec spectral cubes were made spanning the full ob- 
served velocity range. 

In the cross correlation data from runl (confirmed then by 
the analysis of the data from run2) emission was detected at the 
velocities of the known spectral features but also unexpectedly 
in some other spectral channels. The latter was large non noise- 
like side-lobe emission, indicating that maser emission at these 
velocities came from sources located outside of the map region. 
Maps were then iteratively re-centered until the new features, 
were brought close to the image centre. 

The large position offsets of the new features from the cor- 
related position cause large variations of the visibility phase 
over the integration time, and hence introduce amplitude errors 
on the visibilities. These do not affect feature position determi- 
nation, but prevent us from making reliable images. Features F 
& G showed maximal phase changes of 0.6 and 1.8 rarf over the 
integration time of 7 s, producing mean amplitude losses of ~1 
and 8% respectively, on the baseline in Figure[2] 

The positions of the masers were compared to those of 
1 .2 mm dust peak emission that were obtained in a MAMBO- 
2/IRAM-30m map. The image was ob tained in January 2004 
using the MAMBO-2 bolometer array (iRrevsa et alJ ll998) in- 
stalled at the IRAM 30 m telescope in the on-the-fly mode. 
The resulting angular resolution is HPBW ~ 11" and the abso- 
lute pointing is accurate to within ~ 3". 



1 The EVN is a joint facility of European, Chinese, South African 
and others radio astronomy institutes 

2 The Multi-Element Radio Linked Interferometer Network is a 
National Facility run by the University of Manchester, UK 



Fig. 1. General view of the NGC 7538 region. Left: The under- 
lying colour image is a RGB image of J, H, K NIR bands made 
from the 2MASS archive. The grey contours (2, 5, 10, 15, 20, 
30, 40, 50 ,60, 70 ,80, 90 % of peak emission, ~4 Jy beam 1 ) 
show 1 .2 mm dust continuum mapped with MAMBO-2/IRAM- 
30m. White circles indicate the location of the identified in- 
frared sources. White crosses indicate the 6.7 GHz CH3OH 
maser positions. Right: The underlying image is the 1.2 mm 
continuum map. 1 arcmin is 0.78 pc. 



3. Results 

3.1. CH 3 OH masers features F & G 

In the following, we denote the two new spectral features found 
during data reduction as F and G 3 . Their positions and other 
properties are summarised in Tabled 

The comparison of the cross correlated spectrum on the 
shortest baseline with the integrated spectrum of the maser 
emission shows discrepancy in the peak flux of feature F & 
G of ~20% (Pig. Unfortunately, the amplitude errors men- 
tioned in Section [2] and the unreliable baseline subtraction in 
the integrated spectrum introduce errors in the amplitude that 
can be estimated to 15%. These prevent us to invoque resolu- 
tion as responsible for the missing flux. Our best estimate for 
the size of F & G (upper limit) is about half the fringe spacing 
of the shortest baseline, i.e. ~85 mas, or 260 AU, which is to be 
taken as an upper limit. Note that in the higher spectral resolu- 
tion data of run2 (not shown here), F & G show clear multiplic- 
ity (at least two and three components for F& G respectively). 
Observations pointed at the new masers are necessary to obtain 
reliable and accurate maps of all features. 

A review of archival EVN VLBI autocorrelation data for all 
antennas except Effelsberg shows that both F and G are clearly 
present in all experiments since 1997. These features always 
have similar flux ratios to spectral feature A(within a variation 
of the latter of some 10%). This indicates that the emission 
from F and G is probably non-variable. The new features are 
not detected in the autocorrelations from the Effelsberg 100 m 
antenna because they fall outside the HPBW of the primary 
beam at 5cm (i.e. 1.6 arcmin). Note that despite the reduction in 
amplitude due to the primary beam attenuation a hint of feature 
G was present in the discovery spectru m for NGC 75 38 taken 
with the 140 foot Green Bank telescope llMentenll99lh . 

3.2. CH 3 OH masers in NGC 7538 IRS 9 & S 

Figure [0 shows a general view of the NGC 7538 region. The 
CH3OH masers are clearly associated with the most embedded 
sites of high-mass star formation: IRS 1, IRS 9 and S. Masers 
in IRS 1 coincide with the 1.2 mm continuum peak emission 
(within the position accuracy) tracing cold dust emission from 
massive and deeply embedded protoclusters. 

3 The nomencla ture of all features follows the one adopted in 
iMinier et alJ l2000l) . based on the spatial distribution of the masers 
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Fig. 2. Cross correlated spectrum on the KN-DA baseline av- 
eraged over the whole observing run (dashed line) and auto- 
correlation spectrum from the DA antenna (solid line), both 
from runl. Minimal fringe spacing on the KN-DA baseline is 
~190 mas. Features A,B,C,Dand Eare located within 0.5 arcsec 
of IRS 1 . The inset is a zoom-in of the velocity region includ- 
ing features F & G. The multiplicity of F & G is discussed in 
the text. 



Maser feature F coincides with the 1.2 mm dust contin- 
uum emission peak within ~ 2" (Fig. This is less than 
the position accuracy of the 1.2 mm data (~ 3"), and there- 
fore we can state that they are cospatial. The methanol maser 
also coincides with other signposts of high-mass star formation 

in NGC 7538 S as reported in the literature, within the posi- 
■ 1 III I II I 

tion accuracies (Sandell & Sievers 2004; Sandell et al. 2003; 

iKameva etaD Il990ft . In ISandell et alJ J2003I) this protoclus- 

ter/cloud is stated to be one of the most massive of the region, 

with 1000 M Q estimated within a circle of 20"(or 50000 AU) 

in diameter. The authors of that paper also report of a large 

rotating torus of some 14000 AU in radius. Single grey-body 

Spectral Energy Distribution (SED) fitting of the cold dust 

emission in the FIR/(sub)mm/mm emission with an aperture 

of 60" of this source gives a dust temperature of ~35 K and 

an emissivity index (3 — 1.6. This source lies some 10" to the 

south of IRS 1 1 , where a CO outflow is in ferred by the distri - 

bution of the magnetic field orientation jKameva et alJll99ll) . 

On the basis of the above multiple evidence, we consider the 

position of the methanol maser to be the most reliable for the 

protostellar object in NGC 7538 S. 

Maser feature G is coincident within 1" with the peak 

of the 1.2 mm dust emission. The position of G also agrees, 

within the position inaccuracies, with the high resolution 

positions of the peak emission of the continuum at 4.8 

to 10 7 GHz Jvan der Tak et alJ l2000t Ivan der Tak & Mentenl 

2005). Furthermore, one water vapour maser spot as well as 

one class I methanol maser spot seem also to coincide, w ithin 

the positional accu racies, with feature G ( Kam eva et alJ ll990: 

San dell et all2 005). although with an offset in velocity of some 

6kms~'. Because of the fact that class I m ethanol masers are 

mainly found in outflows llKurtz et al.l 2004). we argue that the 



newly detected maser G marks the driving source of the outflow 
(together with the continuum), while the other masers arise in 
th e outflow shocks. Thi s is supported by the geometry proposed 
in Sand ell et alJ (2005) for that source: we might observe the 
outflow from IRS 9 pole-on. The protostellar object in IRS 9 
is accurately located both with the class II methanol maser and 
the cm/mm continuum. 

4. Discussion 

The presence of three methanol maser sources associated to 
three different objects in the same star formation region of- 
fers a unique opportunity to study the evolutionary details of 
methanol maser bearing sources. 

The methanol maser phase is considered to be relatively 
short as compared to the lifetime of the massive protostar, as 
short as a few 10 4 yr (I van der Wall2005T) . This implies that the 
sources hosting methanol masers at 6.7 GHz should show very 
similar properties. The methanol masers in NGC 7538 (features 
A, G & F) on the other hand seem to mark massive young stellar 
objects at different stages of evolution: an IR-quiet protostellar 
object (NGC 7538 S), an IR-bright protostellar object (IRS 9) 
and an UCHu region (IRS 1). The same distinction seems to 
be shown by the flux densities of the masers: 300, 15, 5 Jy for 
A, G, F respectively. 

We propose 2 scenarios for this apparent sequence. The 
first and most intuitive is the evolutionary sequence. Feature 
F might mark the youngest object in the region, feature A the 
most evolved but still within the protostellar phase. This sce- 
nario is also supported by the presence of discs/outflows in the 
maser sourc es: in F a large scale rota ting torus is detected but 
no outflow (Sande lT& Sieversll2004l) : in G a powerful, p rob- 
ably very young, outflow is mapped (Sandell et al. 2005); m 
A a clear disc/outflow geometry is visible (e.g. l Gaume etalJ 
ll995HPestalozzi et aljEooUbe Buizer & Minierl2005l) . More 
evidence has to be collected on the capacity of methanol masers 
to provide such accuracy in the age the sources hosting them. 

The second scenario invokes geometry. The difference in 
1.2 mm dust emission peak flux among the methanol maser 
bearing sources could be due to different viewing angles to the 
sour ce. IRS 1 is seen edge-on, showing the th ick dust in its 
disc (Gaume et al. 1995; Pestalozzi et al. 2004). IRS 9 is sug- 
gested to be seen close to pole-on ( Sandell et al. 2005): the disc 
would be seen almost face-on, not contributing significantly to 
the 1.2 mm flux. NGC 7538 S is still deeply embedded, so dust 
is close to uniformly distributed around the source. This in- 
terpretation seems to be more difficult to accept, because e.g. 
IRS 1 suffers of significant contamination by the very close 
IRS 2 and IRS 3. 

5. Summary and Conclusions 

From the high spatial resolution MERLIN observations of 
6.7 GHz methanol maser emission toward the massive star 
forming region NGC 7538 we have detected and positioned two 
new spectral features. We note that this achievement was possi- 
ble only thanks to the intermediate resolution of the MERLIN 
with its shortest baselines. We conclude the following: 



Pestalozzi et al.: Discovery of two new methanol masers in NGC 7538 



5 



Table 1. Position, line of sight (LOS) velocity (v) and offset (r) 
from feature A of the two maser components F and G. The un- 
certainties in the positions and in the velocity widths are 50 mas 
(135 AU) and 0.044 km s" 1 respectively. 



Masers 


RA 


DEC 


v / FWHM 


r 




J2000 


J2000 


km s -1 


arcmin / pc 


A 


23 13 45.363 


61 28 10.55 


-56.2/1.0 




F 


23 13 44.86 


61 26 48.1 


-53.2/0.5 


1.36/1.06 


G 


23 14 01.78 


61 27 19.7 


-48.8 / 0.9 


2.14/1.67 
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Walsh, A., Burton, M., Hylard, A., & Robinson, G. 1998, 
MNRAS, 301,640 



- on the basis of the accepted idea that methanol masers are 
the most efficient tracers of young and embedded mas- 
sive protostars, the accurate location (to within 50 mas or 
135 AU) of F and G is equivalent with the accurate lo- 
cation of the massive protostars in NGC 7538 IRS 9 and 
NGC 7538 S; 

- the variety of objects marked by methanol maser emission 
could be due to different ages or viewing angle geome- 
try. To test these two scenarios, more accurate observations 
(targeted toward F & G) are required. 
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